Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 5 July 2012 



(MN I^TeK style file v2.2) 



Mind the Gap: Tightening the Mass- Richness Relation 
with Magnitude Gaps 

Andrew P. Hearin 1 , Andrew R. Zentner 1 , Jeffrey A. Newman 1 , Andreas A. Berlind 2 

T ^ , 1 Department of Physics and Astronomy & 

- Pittsburgh Particle physics, Astrophysics and Cosmology Center (PITT PACC), 
University of Pittsburgh, Pittsburgh, PA 15260; 
aphl5@pitt.edu, zentner@pitt.edu, janewman@pitt.edu 

2 Department of Physics and Astronomy, Vanderbilt University, Nashville, TN; a.berlind@vanderbilt.edu 



5 July 2012 



o 

u 

Oh 
6 



> 
o 

o 



X 



ABSTRACT 

We investigate the potential to improve optical tracers of cluster mass by exploit- 
ing measurements of the magnitude gap, m 12 , defined as the difference between the 
r-band absolute magnitude of the two brightest cluster members. We find that in a 
mock sample of galaxy groups and clusters constructed from the Bolshoi simulation, 
the scatter about the mass-richness relation decreases by ~ 15 — 20% when magnitude 
gap information is included. A similar trend is evident in a volume-limited, spectro- 
scopic sample of galaxy groups observed in the Sloan Digital Sky Survey (SDSS). We 
find that SDSS groups with small magnitude gaps are richer than large-gap groups 
at fixed values of the one-dimensional velocity dispersion among group members o~ v , 
which we use as a mass proxy. We demonstrate explicitly that mi2 contains informa- 
tion about cluster mass that supplements the information provided by group richness 
and the luminosity of the brightest cluster galaxy, Lbcg- in so doing, we show that 
the luminosities of the members of a group with richness N arc inconsistent with the 
distribution of luminosities that results from N random draws from the global galaxy 
luminosity function. As the cosmological constraining power of galaxy clusters is lim- 
ited by the precision in cluster mass determination, our findings suggest a new way to 
improve the cosmological constraints derived from galaxy clusters. 
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1 INTRODUCTION 

Galaxy clusters have long been exploited to probe the com- 
position of the universe. The utility of galaxy clusters as cos- 
mological probes using a broad range of tec hniques has been 
reviewed extensively by lAllen et al.l (|201ll ) . Galaxy cluster 
observations are a key component of any effort to constrain 
the cause of cosmolo gical acceleration (see the review by 
IWeinberg et al ] |2012h . Among many advances in cluster cos- 
mology, modern optical surveys have enabled the construc- 
tion of large samples of optically-identified clusters which, in 
turn, have led to competitive cosmological c onstraints from 
optically-identified clu ster abundances (e.g.. lGladders et al.l 
l2007l : lRozo et al.ll201(l ). 

Much of the constraining power of clusters results from 
determinations of their abundance as a function of their 
mass. The abu ndance of clusters by mass may be reliably 
predicted (e.g.. iTinker et al] [2008l. but cluster mass is not 
directly observable. Optical cluster cosmology efforts gen- 
erally rely on using the number of cluster members (rich- 



ness) as a proxy for mass to concurrently fit for cosmologi- 
cal parameters and th e mass-richness relation (although see 
I Newman et al.l (2002) for another technique). In detail, rich- 
ness must be defined precisely for the observational sample 
under consideration, so that the specific definitions of rich- 
ness vary depending upon survey characteristics and cluster 
identification methods. 

One way to improve cosmological constraints from 
optically-identified clusters is to reduce the scatter between 
the observable (richness) and the pre dicted quantity (mass) 
ijRozo et alJl20ld : lRvkoff et al.ll2012l ). In this paper, we sug- 
gest that the differences in absolute magnitude between 
the most luminous cluster members, data already contained 
within optical surveys aiming to perform cluster cosmology, 
can be harnessed to reduce the scatter in cluster mass for a 
fixed set of observables. Specifically, we show that magnitude 
gap, the difference in r-band absolute magnitude between the 
brightest and second brightest members of a galaxy group, 
can aid in the determination of group mass at both fixed 
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richness and fixed r-band luminosity of the brightest group 
member. We provide theoretical and observational support 
for this suggestion usi ng the Bolshoi simu lation of cosmolog- 
ical structure growth (|Klvpin et aJj[201lJ ) an d galaxy group 
and cluster data fro m Data Release 7 of th e Sloan Digital 
Sky Survey (SDSS) (|Abazaiian et al.ll2009l ). In § we de- 
scribe mock galaxy catalogs constructed from the Bolshoi 
simulation. A brief description of the SDSS group data is 
given in § |3] We present results from our mock galaxy cata- 
log in §[4] and from our analysis of the SDSS groups in § [5] 
We draw brief conclusions from these results in § [6] 



scribe a variety of astronom i cal data (see Klypin et al.ll201ll ; 
iTruiillo-Gomez et all l201ll ; IWatson et all |2012| . and refer- 
ences therein). 

Once brightnesses have been assigned to halos accord- 
ing to Eq. ([T]), we construct mock galaxy samples by im- 
posing a brightness cut M r < — 18 on all the mock galaxies 
in Bolshoi. This brightness cut corresponds to Vmax > 92 
km/s . well above the 5 km/s completeness limit of Bol- 
shoi l| Klypin et al.|[201ll ). In this mock catalog, we consider 
galaxy groups to be collections of subhalos associated with 
the same host halo. We study the properties of these mock 
groups in § U 



SIMULATIONS & MOCK GALAXY 
CATALOGS 



We use the Bolshoi N-body simulation |Klvpin et al.|[201ll ) 
to study the connection between magnitude gaps within 
galaxy groups and the mass-richness relation. Bolshoi mod- 
els the growth of structure in a cubic volume 250/i _1 Mpc 
on a side within a ACDM cosmology with total matter 
density Qm = 0.27, Hubble constant h — 0.7, power 
spectrum tilt n s = 0.95, and power spectrum normal- 
ization as — 0.82. The Bolshoi data are available at 
http : //www . mult idark . org and we refer the reader to 



Riebe et alj (|201 lh for additional information. Our anal- 



ysis requires reliable identification of self-bound subhalos 
within virial radii of d istinct halos. We utilize the R0CKSTAR 
ijBehroozi et al.ll201ll ) halo finder in order to identify halos 
and subhalos within Bolshoi. 

To connect the properties of galaxies to dark matter ha- 
los we employ the widel y used subhalo abundance matching 
techn ique (SHAM) (e.g jKravtsov et al.|[2004l ; IConrov et alj 
2006). We assume a monotonic relationship between the 
r-band luminosities of galaxies and the maximum circular 
speeds of test particles within their host dark matter halos, 
Vmax = max GAI(< r)/rj , where M(< r) is the mass of 
the halo interior to the radial coordinate r. In broad terms, 
the reasoning behind this algorithm is that r-band luminos- 
ity is a rough proxy for stellar mass and stellar mass should 
correlate with the depth of the gravitational potential well, 
described by V ma x. Subhalos evolve significantly due to in- 
teractions upon incorporation into a larger distinct halo, so 
Vmax for subhalos may be a poor prox y for stellar mass or 
r-band luminosity. IConrov et alj (|2006h demonstrated that a 
SHAM algorithm that assigns luminosities to subhalos based 
upon the maximum circular speed of the subhalo at the time 
it merged with the distinct halo, V2™ x , can describe a broad 
range of galaxy clustering data from z ~ to z ~ 4. We 
define the circular speed used in our luminosity assignment 
to Vh = Vmax for distinct halos and Vl = Ki^x for subhalos. 
We assign r-band luminosities to halos and subhalos through 
the implicit relation 



n g (> L) = n h (> VL), 



(1) 



where n 9 (> L) is the number d ensity of observed g alax- 
ies with r-band luminosity > L (Bla nton et alj I2005T ) and 
rih(> Vl) is the number density of dark matter halos and 
subhalos with circular speed > Vl- Eq. |TJ ensures that 
the distribution of luminosities assigned to dark matter 
halos and subhalos matches the observed luminosity func- 
tion of galaxies. SHAM models of this kind successfully de- 



3 OBSERVATIONAL DATA 

To study the mass-richness relation observed in low-redshift 
groups and clusters we use a volume-limited catalog of 
galaxy groups identified in Dat a Release 7 of t he SD SS us- 
ing the algorithm described in Berlind et all ((2006) . This 
is an update of the iBerlind et al. I (|2006h groups (based on 
SDSS Data Release 3) to SDSS Data Release 7. All of the 
members of this sample are members of the main galaxy 
sample of SDSS Data Release 7. The group catalog is con- 
structed using a redshift-space friends-of-friends algorithm 
that has been corrected for incompleteness due to fiber col- 
lisions. The particular group catalog we use is constructed 
from galaxies in a volume-limited spectroscopic sample in 
the redshift range 0.02 < z < 0.068 with r-band absolute 
magnitude M r — 51og/i < —19. We refer to this catalog as 
the "Mrl9" group catalog. Each of the 6439 groups in the 
M rl9 catalog contains N > 2 members. We refer the reader 
to lBerlind et al.l (|2006h for further details on the group find- 
ing algorithm. 



4 PREDICTIONS FOR THE POTENTIAL 
UTILITY OF GAP INFORMATION 

In this section we use the mock galaxy catalog described in 
§ [2] to study the predictions of abundance matching for the 
dependence of host halo mass on galaxy group richness and 
magnitude gap. Such a study is idealized for several reasons, 
including the relative simplicity of the SHAM algorithm and 
the fact that groups can be unambiguously identified with 
halos of a particular mass. Nevertheless, this demonstration 
has the distinct advantage that the masses of the host ha- 
los in the Bolshoi simulation are known and do not need to 
be inferred imperfectly from observational data. To proceed, 
we define richness N as the number of mock galaxies within 
the host halo brighter than our specified absolute magni- 
tude threshold. We take the magnitude gap to be the differ- 
ence between the r-band absolute magnitude of the bright- 
est galaxy within the halo and the second brightest galaxy 
within the halo, mia = M r ,i — M r ,2, where M T ,i is the r-band 
absolute magnitude of the i th brightest galaxy in the halo. 

We illustrate our motivation for exploring the utility 
of gap information in Figure [1] in which we plot the mass 
distributions of host halos in Bolshoi in a narrow range of 
richness, 12 < N < 18. The gray, hatched histogram traces 
the mass distribution for all host halos in this richness range, 
while the thin red (thick blue) curve pertains to host systems 
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Figure 1. Plot of the mass distribution of host halos in Bol- 
shoi in the richness range 12 < TV < 18. The mass distribution 
for all host halos in this richness range appears as the hatched, 
gray histogram. The thin red curve traces this distribution for 
the subsample of host systems with a very large magnitude gap 
(mi2 > 1.5). The mass distribution of host halos with m\2 < 0.2 
is plotted with the thick blue curve. Note that most objects are 
not in either m\2 bin traced by the red or blue histograms, which 
is why the full distribution plotted in gray does not resemble a 
combination of the red and blue. Evidently, large gap systems 
tend to be more massive at fixed richness, suggesting the pos- 
sibility that mi2 can be used to improve the calibration of the 
mass-richness relation. 



with a large (small) magnitude gap mi2. Evidently, large 
gap systems tend to be more massive than their small gap 
counterparts at fixed richness. Of course the trend in Fig. [1] 
could simply be due to the finite width of the richness bin we 
have chosen, rather than following as a consequence of the 
mass-richness residual being correlated with mi2. To explore 
this issue more rigorously, we employ standard regression 
analysis techniques to find the linear relationship between 
ln(jV) and ln(M) that minimizes cr(ln M). The use of a linear 
regression is well-m otivated by previous results (for example, 
iBecker et all 120071 ) that find the mass-richness relation to 
be well-described by a power law. We find that our best 
fit model, M fit (iV) = CN a , with C = 2.2 x 1O 12 M , and 
q = 1.1, gives an accurate description of the mass-richness 
relation for the rich groups (N > 10) in our mock sample, 
yielding a mean residual (5 InM) ~ 0.02, and a residual 
dispersion of <x(lnM) ~ 0.33. 

In Figure [2] we plot the mean residual 8 In M as a func- 
tion of mi2. The trend suggested by Fig. [T] is borne out: 
groups and clusters with a large (small) magnitude gap mi2 
are more (less) massive than the average M at a given rich- 
ness. A linear fit to the results illustrated in Fig. [2] indicates 
that 8 InM oc 0.18mi2, implying that there may be signif- 
icant information about a cluster's mass contained in the 
magnitude gap mvz that is not contained in the richness 
alone. Moreover, the mass-dispersion about the mvz — N 
plane determined by a two-dimensional linear regression im- 
proves by 18% to aQxiM) ~ 0.27, further demonstrating the 
potential improvement in mass determination provided by 
the use of magnitude gap information. 

While the trend in the mass-richness relation with mi2 



Figure 2. Plot of the residual 8 In M of the best-fit mass-richness 
power law relation exhibited by groups and clusters in Bolshoi as 
a function of mi2. A linear fit to the residuals indicates that 
8 InM oc 0.18mi2, suggesting that exploiting the magnitude gap 
may significantly improve cluster mass estimation techniques that 
rely only on richness. 



as seen in Fig. [2] is the novel feature of this work, the ten- 
dency for large gap systems to be more massive at fixed 
richness appears in a variety of guises in the literature on 
fossil groups, which are typically defined to be galaxy groups 
with mi2 > 2. The picture of fossil groups that is the most 
prevalent in the literature is that these are groups which 
assembled most of their mass at high redshift, so that pro- 
cesses such as dynamical friction and mass loss have had 
ample time to deplete these systems of their most mas- 
sive satellite galaxies, leaving behind a very b right central 
galaxy with few comparably bright satellites (jjones et al.l 
l2003l : jp'Onghia et al.ll2005l : (zentner et aD l2005| ) . 

In the above scenario for the origins of fossil groups, the 
same processes that lead to the formation of a large magni- 
tude gap are also at work in the depletion of the number of 
group members above a given brightness threshold. In other 
words, this picture posits a dynamical connection between 
mi2 and richness. However, even in the complete absence of 
such dynamical processes, we may still expect large-gap sys- 
tems to have fewer members than small-gap systems. This 
is a consequence of the shape of the Schechter function: for 
any luminosity function <&(!/) with a slope that steepens 
with brightness, the average gap 77112 obtained from a set of 
TV random draws from $(!/) increases as N decreases. 

The extent to which dynamical processes influence the 
relationship between 77J12 and N remains an open question, 
but we point out here that if 77112 were purely statistical, 
resulting from TV random draws from a Schechter function, 
the magnitude gap would contain no information about mass 
that would not already be provided by knowledge of richness. 
This follows from Bayes' Theorem. If mn were determined 
strictly by N random draws from a common luminosity func- 
tion, then P(m 12 \M,N) = P(m 12 \N). Then it follows di- 
rectly from Bayes' Theorem that P(M\N,m 12 ) = P(M\N); 
that is, the probability distribution of mass is unchanged by 
knowledge of 77112 when the richness is known. This is sen- 
sible because in this scenario the 77112 distribution is deter- 
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mined entirely by N, so knowledge of m\i provides no new, 
independent information about the system. Thus the trend 
exhibited by the mass-richness relation seen in Fig.[2]reflects 
a relationship between the magnitude gap and the mass of 
host halos in our mock sample beyond that expected from 
sampling a luminosity function a finite number of times. 

The SHAM-based results presented in this section show 
that a simple and well-supported model for populating ha- 
los with galaxies makes a definite prediction for the rela- 
tionship between halo mass, magnitude gap, and richness in 
groups. To be sure, there are many factors that will tend 
to wash out the clear trend seen in Fig. For example, 
group/cluster membership (and hence, richness) is deter- 
mined by halo membership, so contamination by interlopers 
due to projection effects is not included; the simulations do 
not suffer from edge effects that may alter the richnesses 
of groups in real surveys; and the fidelity with which the 
SHAM prescription for identifying halos and subhalos with 
luminous galaxies correctly describes richness has not been 
extensively tested in this context. As we will see in the next 
section, despite these complicating factors the dependence 
of the mass-richness relation on mi2 appears to be signifi- 
cant in a spectroscopic sample of galaxy groups observed in 
SDSS. 



5 THE OBSERVED CONNECTION BETWEEN 
MASS, RICHNESS, AND MAGNITUDE GAP 

We now examine the mass-richness scaling relation seen in 
groups and clusters in the volume-limited Mrl9 galaxy group 
catalog described in § [3] The principal result of this section 
appears in the bottom right panel of Figure [3] which we 
will argue demonstrates that the magnitude gap contains 
information about cluster mass that is independent from 
both richness and Lbcg- Most of the work described in this 
section, illustrated in the remaining three panels of Fig. [3] 
addresses a variety of possible systematics and selection ef- 
fects that are germane to our primary conclusion. 

As our mass proxy for the groups we use a v , the one- 
dimensional velocity dispersion of member galaxies, defined 
as 



1 + 2 



N 



1 ™ 



(2) 



where N is the number of group members, Zi are the red- 
shifts of the member galaxies, and z is the redshift of the 
(unweighted) group centroid. We are chiefly interested in de- 
termining whether the scaling of richness with a v changes 
when comparing samples of groups with different mi2, as is 
suggested by the results in §|4] Assigning mi2 values to the 
groups requires some care due to complications presented by 
fiber collisions which affect ~ 8% of the galaxies in Mrl9. 
Fiber collided galaxies in the Mrl9 sample are assigned the 
redshift of their nearest neighbor on the sky. As a galaxy's 
redshift is used to infer its absolute magnitude, this can 
result in a large error in the inferred brightnesses of fiber 
collided galaxies. Thus if either the brightest member of the 
group (BCG) or the second-brightest member (which we de- 
note below as the SBG for 'second-brightest galaxy') is fiber 
collided, the group will be assigned an erroneous value of 



mi2. As the BCG in particular is likely to be found near 
the group centroid where the galaxy number density can be 
very large, this scenario is relatively common; we find that 
~ 2% of the groups in Mrl9 have either a BCG or a SBG 
that is affected by fiber collisions. We define mi2 to be the 
magnitude difference between the two brightest non-fiber- 
collided galaxies in the group, but we note that all of our 
conclusions remain unchanged if we allow galaxies affected 
by fiber collisions to be used in the definition of the magni- 
tude gap. 

In the upper left panel of Figure [3] we plot the mean 
richness of galaxy groups in Mr 19 as a function of a v . This 
panel illustrates the results of a calculation in which we have 
divided the Mrl9 groups into logarithmically-spaced bins 
over the range 150km/s < a v < 350km/s, and computed 
the mean richness in each bin. Evidently, large-gap systems 
do indeed appear to be less rich at fixed a v (mass) than 
small-gap systems. However, before drawing this conclusion 
we now address several possible selection effects and com- 
plicating factors that could influence this result. 

Selecting groups with a large magnitude gap biases 
one to select groups with a luminous BCG. For example, 
the Mrl9 galaxy sample contains no galaxies dimmer than 
M r = —19, so selecting groups with mi2 > 1.5 requires 
the BCG to have a brightness M r < —20.5. Meanwhile, 
iReves et al.l (|2008l ) suggested that BCG luminosity provides 
information on group mass that is independent of richness, 
so we must account for this bias to ensure that the magni- 
tude gap provides information that is independent from the 
known correlation between mass and BCG luminosity. 

We illustrate the differences in BCG luminosity induced 
by this selection in the lower left panel of Fig. UJ with the 
thin blue (thick red) histogram tracing &bcg(L) for small- 
gap (large-gap) systems. The upper right panel of Fig. [3] 
is an illustration of the potential importance of this effect 
as BCG luminosity clearly informs the a v — N relation of 
the SDSS galaxy groups. Systems with brighter BCGs are 
also richer at fixed a v . This trend has the sense that should 
be expected if Lbcg were determined by the brightest of 
N ra ndom draws from a fixed , global luminosity function 
fe.g.. IParaniape fc Shetr]|201lf ). While there may be other, 
dynamical effects at work in galaxy groups that influence the 
correlation between Lbcg and group mass, such processes 
do not counteract the positive correlation between Lbcg 
and richness at fixed a v seen in our group sample. We will 
return to this issue shortly. 

To account for bias from differences in the BCG lu- 
minosity distribution, we have drawn a random subsample 
of 1000 of the low-gap, mi2 < 0.2, groups with a BCG 
brightness distribution that matches that of the large-gap, 
mi2 > 1.5 groups. The magenta, hatched histogram in the 
bottom left panel of Fig. [3J shows <J?bcg (L) of the result- 
ing LBCG-matched subsample of the m\i < 0.2 groups. The 
a v —N scaling relation for the small-gap, matched-BCG sys- 
tems is shown as the purple triangles in the upper left panel 
of Fig. [3j Eliminating any differences between the BCG lu- 
minosity distributions of low-gap and high-gap systems in- 
creases the disparity between the a v — N relations of low- 
and high-gap groups, so mi2 informs this relation in a man- 
ner that is independent from BCG luminosity alone. 

We compare the a v — N scaling relation for large-gap 
systems and their matched-Z/BCG, small-gap counterparts 
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in the bottom right panel of Fig. [3] The principal result of 
our analysis in this section lies in the comparison between 
the red diamonds and magenta triangles: at fixed Lbcg, 
large-gap systems are under-rich relative to small-gap sys- 
tems at fixed a v (mass). However, because of the natural 
correlation between m\2 and N discussed in § U some care 
is required before interpreting the differences between these 
points as implying that large-gap and small-gap systems 
have intrinsically different mass distributions at fixed rich- 
ness. The difference between the red diamonds and magenta 
triangles shows that P(N\o~ v , 77112) 7^ P(N\a v ), or equiva- 
lently, P(mi2\o v , N) 7^ P(m\2\o v ). This inequality would 
hold even in a universe in which mi2 is determined solely by 
TV random draws from a global luminosity function. Yet, as 
discussed in § [4] in such a universe the magnitude gap con- 
tains no information about cluster mass that is independent 
from richness. However, it is possible to show that the o~ v —N 
scaling relations exhibited by the SDSS groups are distinct 
from the relations expected if mi2 were solely determined 
by random selection from a global luminosity function. 

We demonstrate that this is the case by contrasting the 
SDSS group data against the following Monte Carlo (MC) 
simulation. We draw a random galaxy group from our SDSS 
sample and assign the values of a v and N of this group 
to the "MC" galaxy group. We populate the MC group by 
drawing N galaxies from the global luminosity function of 
all the galaxies in the Mrl9 group sample. We assign the 
MC group a value of m\2 by taking the difference between 
the r-band absolute magnitudes of the brightest and next- 
brightest galaxies used to populate the MC group. We repeat 
this procedure 10 6 times to construct a sample of one million 
MC groups with the same group multiplicity function and 
o~ v — N scaling relation as the Mrl9 sample. The MC groups 
exhibit a correlation between m\2 and N at fixed a v that 
is determined solely by the statistics of random draws from 
the global luminosity function. In the bottom right panel of 
Fig. [3] we plot the a v — N scaling relation for MC groups 
with mi2 > 1.5 with orange asterisks, and MC groups with 
mi2 < 0.2 with black crosses. Both samples of MC groups 
have been selected so that their Lbcg distributions match 
that of the large-gap systems seen in the data, so that all 
the samples plotted in the bottom right panel have the same 
$>bcg(L). 

The differences in richness at fixed o~ v between the large- 
gap and small-gap groups in the SDSS data are clearly more 
significant compared to their MC counterparts. Based on 
this exercise, we can conclude that the gap distribution is 
not determined solely by statistical effects, providing strong 
evidence that having a large magnitude gap is, in fact, cor- 
related with being under-rich at a given mass. Moreover, 
because we have controlled for the brightness of the BCG, 
this implies that the magnitude gap mi2 contains informa- 
tion about the masses of galaxy groups that is independent 
from both richness and BCG brightness. These results sug- 
gest that it may be possible to use mi2 to improve optical 
estimators of group and cluster masses. 

As a further check that the magnitude gap provides 
information that is independent from richness, we have sup- 
plemented the comparison of the Mrl9 groups to the MC 
groups by comparing the mi2 distribution of two subsam- 
ples of the Mrl9 groups: one with 75 km/s < o~ v < 150 km/s, 
and the other with 200 km/s < a v < 400 km/s, both with 



the same distribution of richnesses. We find that at fixed 
richness, the groups with larger velocity dispersion have a 
higher proportion of large mi2 values relative to the small 
a v groups such that the mean 77112 values of the two dis- 
tributions are offset by ~ 2.5<r. The sense of this trend is 
in keeping with the results presented in Fig. [3] This result 
explicitly demonstrates that P(mi2\o~ v , N) 7^ P(mi2\N), di- 
rectly implying that P(a v \mi2, N) 7^ P(a v \N), i.e., that the 
magnitude gap is informative about group/cluster mass even 
when the richness is known. 

We conclude this section by addressing an additional 
possible systematic error in our data analysis. In our mock 
sample, the typical scatter in <j v at fixed host mass is ~ 15 
km/s, and so we do not expect bin-to- bin scatter in mass to 
be a significant source of error. However, the two samples 
plotted in the bottom right panel of Fig. [3] have different 
richness distributions, and so their velocity dispersions a v 
are not determined with equivalent accuracy. In particular, 
as the error in o v depends upon N, this may induce system- 
atic differences between the measurements of the velocity 
dispersion of the groups in the two samples, potentially pro- 
ducing a spurious difference in the o v — N relation exhibited 
by small-gap and large-gap systems. To estimate the signifi- 
cance of this effect we have conducted the following exercise. 
For every group in each sample, we randomly select four 
members and use only these members to compute the veloc- 
ity dispersion. We refer to the dispersion thus computed as 
the reduced velocity dispersion, <rjj ed . For the small-gap sys- 
tem we find that the slope of the a v — N relation is slightly 
steeper than the al cd — N relation (as expected, since the 
cr„ measurements are generally noisier than a v estimates), 
but the trend with gap persists at similar levels and our 
conclusions remain unchanged: at fixed o"Jf d , groups with a 
large magnitude gap have fewer members than groups with 
small values of mi2- 



6 SUMMARY & CONCLUSIONS 

We have studied the improvement to group and cluster mass 
determination that may be reaped by exploiting the magni- 
tude gap, mi2, between the two brightest group and clus- 
ter members in addition to group richness, N. After fitting 
the mass-richness relation in our mock sample of simulated 
groups and clusters with a power law, we find a signifi- 
cant correlation between the magnitude gap and group mass 
residual, S InM oc 0.18 m\2- The strength of this correlation 
is significant when compared to the scatter about our power 
law fit, a (5 InM) = 0.33. 

We see a similar trend in a volume-limited spectroscopic 
sample of galaxy groups observed in the SDSS. For group 
samples with different magnitude gaps, we find that large- 
gap grou ps have fewe r mem bers than small-gap groups at 
fixed cr„. iReves et al.l (|2008T ) recently suggested that using 
the luminosity of the BCG could aid in reducing scatter in 
cluster mass determinations. Similarly, we found that BCG 
luminosity is correlated with richness at fixed o v in our SDSS 
groups. By constructing appropriate random samples from 
the SDSS data, we were able to conclude that m\2 contains 
information about group mass that is not contained in either 
richness or Lbcg- 

Our findings are in conflict with recent results suggest- 
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Figure 3. Upper Left: Richness as a function of mass proxy, a v . Points are placed at the median a v in each bin. Horizontal error bars 
represent the inner quartiles of the a v distribution in each bin. Vertical error bars show the error on the mean. Red diamonds show 
systems exhibiting a large magnitude gap m\2 > 1.5, while blue squares show small-gap systems. The magenta triangles represent a 
randomly selected subsample of small-gap systems with a BCG luminosity function ^BCGf^i matched to that of the large-gap systems. 
Lower Left: Thick (Thin) histograms show the distribution of BCG luminosity "J'bcgC^O (normalized to unit area) of large- (small-) gap 
systems in Mrl9. The magenta, hatched histogram traces the BCG luminosity function of a random subsample of small-gap systems 
selected to match $bcg (•£) °f the large-gap systems. Upper Right: Mean richness as a function of a v for galaxy groups with different 
BCG brightnesses. Lower Right: The cr v — N scaling relations for large-gap groups and small-gap groups. All samples have the same 
&bcg(L)- Orange asterisks and black crosses represent results from Monte Carlo (MC) simulation of the group population in a model 
assuming that gap does not explicitly depend upon a v (the group mass proxy). MC groups have negligible errors on the mean richness. 
Red diamonds and magenta triangles represent large- and small-gap groups in the SDSS data, as in the upper left panel. Differences 
between large-gap and small-gap groups in the data are significant compared to the MC groups, demonstrating that the magnitude gap 
mi2 contains information about group/cluster mass that is independent from both richness and I/bcg- 



ing that the distri bution of magnitude gaps is determined 
purely by richness. IParaniape fc Shethl f|201 lh find that the 
abundance of galaxy groups in as a function of mag- 

nitude gap is consistent with the distribution resulting from 
a set of random draws from a global luminosity function, 
implying that group mass is only related to m\i through 
mutua l covariance with richness. Meanwhile, IProctor et al.l 
(|201lf) claim that the fossil fraction, defined as the frac- 
tion of groups with ran > 2, is purely a reflection of the 
abundance of low-richness systems. In demonstrating that 
groups with different magnitude gaps exhibit different re- 
lationships between a v and richness we have established 



1 Note th at the effective volume of th e Mr 19 galaxy group sample 
studied in IParaniape fc Shethl j201lT) . which was based on SDSS 
Data Release 3, is less than half the effective volume of our Data 
Release 7-based Mrl9 sample. 



that P{m\i\a v , N) ^ P(m\2\N), explicitly showing that the 
magnitude gap is not strictly determined by richness. 

We have demonstrated that, in principle, magnitude 
gap can inform cluster mass determination both in simple 
mock catalogs constructed from N-body simulations and in 
spectroscopic SDSS data; however, we leave the determina- 
tion of the extent to which these relations may aid forthcom- 
ing cluster cosmology efforts as a subject of future work. 
Existing cluster cosmology samples differ from the group 
catalogs we hav e studied in several ways. For example, the 
maxBCG clusters l|Koester et al.ll2007l ') are selected using pho- 
tometric (rather than spectroscopic) data with an indepen- 
dent algorithm that uses color and i-band luminosity crite- 
ria. Additionally, the maxBCG clusters extend to higher lu- 
minosities and larger richnesses than our groups. Forthcom- 
ing cluster cosmology efforts with data from imaging surveys 
like the Dark Energy Survey (DES) will likewise identify 
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clusters using photometric data and probe larger richnesses. 
Moreover, the richnesses of maxBCG clusters are defined 
according to a more compl ex optimization procedure than 
the r ichnesses of our groups dKoester et al.ll2007l:lRc>zo et al.l 
l2009h . We have analyzed the iBerlind et al.l (|2006l ) clusters 
because the cluster membership assignments for the maxBCG 
clusters are not publicly available. We do not anticipate the 
correlations that we point out here to be particularly strong 
functions of redshift or richness, but this will need to be 
tested more extensively both in mock catalogs and forth- 
coming data. 

There are at least two distinct ways in which m\i 
may be exploited to tighten the relationship between 
group/cluster mass and richness. First, the magnitude gap 
could be treated on an equal footing with richness, so that 
rather than calibrating the mass as a one-dimensional func- 
tion of richness one could instead treat the mass as a func- 
tion of N and mi2 simultaneously. Simulations coupled with 
detailed studies of extant and near- future data could provide 
parameterized forms for the mi2 — M relation with reason- 
able priors as they do now for the N — M relation. We stud- 
ied the potential benefit of this approach in our mock group 
sample by comparing the difference in the scatter about the 
residual mass estimation between a one-dimensional linear 
regression on richness and a two-dimensional linear regres- 
sion on richness and magnitude gap. We find that the scatter 
in the residuals 6 lnM improves by ~ 15 — 20% when us- 
ing a fit for mass as a function of N and 777,12 instead of N 
alone. While this improvement may seem modest, it comes 
at no additional observational cost, because the magnitude 
gap will always be available in the same data set used to 
measure the richness. 

A second approach is suggested by the observation that 
the relationship between 5 lnM and mi2 in our mocks ap- 
pears to be nonlinear, with the systems with the very largest 
gaps appearing to be outliers in the mass-richness relation 
(see Fig. [2]). The highest-gap systems have inordinately large 
masses at fixed richness. One may imagine utilizing gap in- 
formation to identify significant outliers in the mass-richness 
relation. It may be possible to impose a cut on mi2 and 
restrict consideration to groups and clusters with a modest 
magnitude gap (mi2 < 1.5) to calibrate the mass-richness re- 
lation. This may be particularly helpful in photometrically- 
identified groups because interloper contamination will be 
more significant in the absence of spectroscopic redshifts, 
but interlopers can only reduce the magnitude gap, so large- 
gap systems will remain mass-richness outliers. Of course, 
the effect of any such cut on cluster abundance must be cal- 
ibrated and accounted for, and this must be a subject for 
further work. In a forthcoming companion paper (Hearin et 
al. 2012, in prep), we study in detail the global abundance 
of groups and clusters as a function of mi2 over a wide range 
of masses, providing precisely the information that would be 
required to carry out this second approach. 

As we were preparing to submit this paper for publi - 
cation, we became aware of another paper l|Wu et al.ll2012l ) 
that studies some of the same material that we have ad- 
dressed. In particular, they focused on the host halo with 
the most extreme difference between the Kn ax value of the 
host and its largest subhalo, the N-body simulation analog 
of a cluster with a very large magnitude gap. They found 
that this halo also appeared to be an outlier in many of 



the host halo properties they studied, including the number 
of subhalos contained by the host. This finding appears to 
be in keeping with the second approach described above to 
using mi2 in the calibration of the mass-richness relation. 

However the calibration is conducted, our results sug- 
gest that the magnitude gap 777,12 contains information about 
cluster mass that is independent from both richness and 
Lbcg, the luminosity of the brightest cluster member. Ex- 
ploiting this additional information to improve existing op- 
tical tracers of cluster mass may improve the constraining 
power of optically-identified galaxy clusters on cosmology. 
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